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Existence of a protected surface state described by a massless Dirac equation is a defining property 
of the topological insulator. Though this statement can be explicitly verified on an idealized flat 
surface, it remains to be addressed to what extent it could be general. On a curved surface, the 
surface Dirac equation is modified by the spin connection terms. Here, in the light of the differential 
geometry, we give a general framework for constructing the surface Dirac equation starting from the 
Hamiltonian for bulk topological insulators. The obtained unified description clarifies the physical 
meaning of the spin connection. 



I. INTRODUCTION 

The emergence of low-energy electron states on the sur- 
face of a three-dimensional (3D) strong topological insu- 
lator is topologically protected by a nontrivial value of 
the strong Z2 index (i.e., Z2 = These states on a 

flat surface obey a massless Dirac equation, and possess a 
linear energy dispersion forming a gapless cone-like struc- 
ture (i.e., Dirac cone) in reciprocal space. The existence 
of a single Dirac cone has been verified experimentally by 
a number of angle-resolved photoemission spectroscopy 
measurements performed in different realizations of the 
3D topological insulator 4-8 and is now incontrovertible. 9 
It should be noted that these features of Dirac electrons 
manifest themselves on the flat surface of topological in- 
sulator samples. There arises a natural question: how do 
low-energy electrons behave on a curved surface? The 
closed surface of any realistic topological insulator sam- 
ple is inevitably curved, and we expect that this curved 
nature plays a role in some cases. 

The behavior of low-energy electrons on a curved sur- 
face of, such as, cylindrical and spherical topological in- 
sulators 1 ^— is of considerable theoretical interest. The 
case of a cylindrical topological insulator— ~ 14 i 16 i 19 has 
been studied by the use of an appropriate Dirac equation. 
A remarkable feature is the appearance of a finite-size en- 
ergy gap separating the upper and lower Dirac cones of 
surface electron states. This is the consequence of half- 
integral quantization of angular momentum caused by a 
Berry's phase tt. The origin of the Berry's phase has 
been explained as follows: since the real spin of surface 
electrons is locked on the tangential plane at any point 
on a carved surface (i.e., spin-to-surface locking), surface 
curvature induces a Berry's phase tt on going around the 
cylindrical surface. It has been pointed out that the same 
type of the Berry's phase tt appears in the spherical case, 
and therefore a finite-size energy gap also appearsj 10 ' 15 i 18 

In spite of the accumulation of such studies, we have 
so far not reached the unified theoretical framework for 
analyzing the behavior of the low-energy electrons on 
a curved surface. It is believed that low-energy states 
are inevitably described by a Dirac-type equation even 
on curved surfaces. However, an explicit demonstration 
of this has been lacking. We here briefly summarize 



the present status of the theoretical studies focusing on 
curved surfaces of topological insulators. The analyses on 
this problem have been so far performed mainly for cylin- 
drical and spherical systems, to which the following two 
approaches have been applied. The first group of stud- 
ies is based on the two-dimensional (2D) Dirac equation 
for a flat surface, and takes account of the curved na- 
ture of the surface by a coordinate transformation! 1 ^— 
Certainly, the resulting curved surface Dirac theory can 
be applicable to an arbitrary curved surface in this ap- 
proach, but for a clear reason that it assumes a 2D Dirac 
theory from the outset, it fails to answer the question 
whether the low-energy electrons on an arbitrary curved 
surface obey indeed a Dirac-type equation. Furthermore, 
this has a drawback that it ignores the 3D nature of the 
original problem. In the second category of approaches 
one starts from a 3D bulk Hamiltonian, and an effective 
2D theory for a given curved surface is derived by the use 
of the k ■ p approximation* 1 ^— This can in principle take 
account of the 3D nature. In this paper we generalize the 
second approach and construct a unified framework for 
analyzing the curved surface problem of topological in- 
sulators. 

Construction of the general framework reveals the fol- 
lowing issues related to the curved surface Dirac theory. 
One is an insufficiency of the argument on the Berry's 
phase tt. The Berry's phase tt is recognized as a con- 
sequence of the so-called spin-to-surface locking in liter- 
atures jiLLr— ii& However, since this locking breaks down, 
for example, in spherical topological insulators^ the ar- 
gument based on the spin-to-surface locking is not more 
than a phenomenological interpretation. Indeed, the ap- 
pearance of the Berry's phase tt should be attributed to 
something more fundamental. Another issue to be ad- 
dressed concerns how one should interpret the spin con- 
nection in the curved Dirac theory. An effective Hamil- 
tonian in the curved Dirac theory contains a fictitious 
vector potential that is often referred to as the "spin 
connection" . The spin connection is a mathematically 
well-defined concept, and encodes the information on the 
curved nature of the surface. Yet, its physical origin re- 
mains mysterious, at least for those researchers in the 
condensed matter community. Here, in this paper we 
attempt an accessible presentation of this issue. 
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FIG. 1. Curved surface of a 3D topological insulator: cross 
section in the plane spanned by ei and e^. 



In the next section we present a 3D model Hamilto- 
nian for a strong topological insulator and define a set 
of curvilinear coordinates that can be used to describe 
an arbitrary curved surface. In §3, we introduce a gen- 
eral framework to derive the effective surface Hamilto- 
nian for the low-energy electrons within the k ■ p approx- 
imation. We show that the effective Hamiltonian is in- 
deed expressed in a generalized Dirac form with the spin 
connection. The role of the spin connection is clarified 
and given an intuitive interpretation. In §4, we discuss 
the boundary condition for a spinor wave function in spa- 
tially periodic systems, and identify the precise origin of 
the Berry's phase tt. In §5, our framework is applied to 
two known examples to demonstrate how it works. The 
last section is devoted to summary and discussion. We 
set H — 1 throughout this paper. 



II. MODEL 

Let us start with the following bulk effective Hamilto- 
nian for a 3D isotropic (strong) topological insulator in 
the continuum limit i 21 ' 22 i?buik = Tn(p)r z + At x (<t ■ p), 
where p = — iV and m(p) = mo + "i2P 2 is the mass term. 
Without loss of generality, we assume that mo > and 
m 2 < 0. The two types of Pauli matrices cr — {a Xl a y ,a z ) 
and t = (t x ,t v ,t z ), respectively, represent the real and 
orbital spin degrees of freedom. If the ordinary matrix 
representation of r is used, i?b u ik is expressed as 



bulk 



m(p) A(tr ■ p) 
A(cr ■ p) —m(p) 



(1) 



Let us consider a curved surface of the sample de- 
scribed by two coordinates (x 1 ^ 2 ) as X a — f a {x 1 1 x 2 ) 1 
where X = (X 1 , X 2 , X 3 ) represents the position of an 
arbitrary point on the surface in the 3D Cartesian coor- 
dinates (see Fig. 1). Depending on the geometry (e.g., 
closed vs. open surfaces, etc.), the coordinate x\ repre- 
sents either a linear (non-cyclic) or a cyclic coordinate. 



For simplicity we focus on samples whose entire surface 
is described by a single set of functions {/"}, although 
with little modification our argument can be extended 
to cases where the entire surface is wrapped by several 
patches and a separate set of functions is needed to de- 
scribe each patch. We introduce a curvilinear coordinate 
system which is suitable to analyze low-energy electron 
states of the topological insulator localized near its sur- 
face. It should be noted that as surface electron states 
has a finite penetration depth A, we require that a curvi- 
linear coordinate system is well-defined only in the sur- 
face region of width on the order of A. Let e\ and e 2 be 
the two tangent vectors defined by 



dx l 



(2) 



Note that e\ and e 2 are not necessarily orthogonal with 
each other nor normalized to be unity. Let e 3 be the unit 
normal vector defined by 



e.3 



ei x e 2 
lei x e 2 | 



(3) 



which for simplicity is assumed to be outward normal 
to the surface. We introduce the third (perpendicular) 
coordinate x 3 along the straight line designated by e 3 , 
and set a; 3 = x s { just on the surface. Let us consider fic- 
titious internal surfaces obtained by varying (a; 1 ,^ 2 ) at 
fixed values of x 3 satisfying x s f > x 3 > x s f — A. These 
surfaces are well-defined over the entire sample as long as 
the smallest value of local radius of curvature at x 3 = x s { 
is much longer than A. Let f a {x 1 : x 2 1 x 3 ) be the function 
that describes the internal surface at x 3 . In terms of this 
we define ei and e 2 similar to eq. ©. Now ei and e 2 be- 
come functions of (a; 1 , x 2 , x 3 ), while e% is independent of 
x . We employ ( 3 ) as the curvilinear coordinates 

with the basis vectors 



ei = ei(x ,x ,x ), 
e 2 = e 2 {x 1 ,x 2 ,x 3 ), 
e 3 = e 3 (x 1 ,a; 2 ). 

Let us introduce e 1 that satisfies e,- • e J = 5s 



(4) 
(5) 
(6) 

We 



also introduce the metric tensors in a symmetric bilinear 
form defined by gij = e; • ej and g lJ = e l ■ e? , which 
satisfy gik9 = Here and hereafter we use the con- 
vention that a repeated index, such as i, j, and k, should 
be summed over 1,2,3. Obviously, #33 = g 33 = 1 as 
|e 3 | = 1, and g 13 = ,g 23 = g 13 = g 23 = as ei_Le 3 
and e 2 _Le 3 . In this coordinate system, an infinitesimal 
volume element is given by dV = dSdx 3 with 



dS = VGdx^x 2 
being an infinitesimal area element, where 
^(x 1 ,^ 2 ,^ 3 ) = det{g l3 } = gng 2 2 - 



9l2- 



(7) 



(8) 



Within the curvilinear coordinates presented above, the 
3D Cartesian coordinates (X 1 , X 2 , X 3 ) — (x,y,z) and 
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the differential operators with respect to them are rep- 
resented as X a — eiaX 1 and d/dX a = e l a di, where 
di = d/dx 1 . For later convenience we decompose the 
Laplacian into two parts as V 2 = Aj_ + Am, where 



A, = ±B> (VGd 3 ) 



(9) 
(10) 



Finally we present a set of equations which describe 
the spatial variation of e, against x 1 and x 2 . The tangent 
vectors obey the Gauss equation, 



fe=i 



where 



hj = ~(dje 3 ) ■ e i; 



r 



(11) 

(12) 
(13) 



The unit normal vector obeys the Weingarten equation, 
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a,-e 3 = -5^6/e i) (14) 



where 



V UKC:;:-C. (15) 

fe=l 

These equations are used in Appendices B and C in de- 
riving, or simplifying, matrix elements of the effective 
Hamiltonian. 



III. DERIVATION OF THE EFFECTIVE 
HAMILTONIAN 

Let us now derive the effective surface Hamiltonian 
for low-energy electron states of the topological insulator 
for which curvilinear coordinates (x , x 2 , x 3 ) are defined. 
The derivation consists of four steps. 

In the first step, we rewrite eq. (H} in terms of the curvi- 
linear coordinates, and then divide it into perpendicular 
and parallel components; the former describes the pene- 
tration of surface wave functions into the bulk while the 
latter determines low-energy properties of surface states. 
As a result, -ffbuik is decomposed as ifbuik = H±_ + Hu 
with 



m 



rri2A± -iAa 3 d 3 



-iAa 3 d 3 
-m 2 A|i 



-m + m 2 Aj_ 



m2A|| 



(16) 
(17) 



where we have used er ■ V = a 1 ■ di with 



(18) 



In the second step we solve the perpendicular eigen- 
value equation Hx\tp) — E^\ip) j 16 ! 18 , 21 ! 22 anc j obtain two 
basis states. The appropriate boundary condition for 
is \ip(x 3 = x s [)) = 0. That is, all four components of the 
wave function vanish on the surface at x 3 — x s {. As 
the simplest approximation, we replace A^ in H± with 



A± = d 2 



^(d 3 lnG)d 3 
d 2 + A(x\x 2 )d 3 , 



(19) 



where the definition of the average (• • • ) over x 3 is given 
below [see eq. (1331) ]. Then, we can show that the eigen- 
value equation has surface solutions of the damped form, 
\tp ) = e K ( x ~ Xet '\u). Here k characterizing the penetra- 
tion depth A is determined by det{M^} = with 



Mi 



where £ 



mo — TO2C — E± 
—\Ak<j 3 



-\Akcj 3 
-too + m 2 ( - E\ 



, (20) 



-An. Let us introduce the two eigenvectors 



n± of a . They satisfy 



a 3 n± 



±n±, 



(21) 



and are regarded as local spin quantization axis. It 
should be stressed that n± points in the ±e 3 direction if 
the spin axes (s x , s y ,s z ) are identified with the 3D Carte- 
sian coordinates (x,y,z). That is, n + is outward 
(inward) normal to the tangential plane at (x , x 2 ) on the 
surface. In terms of the 4x4 unitary matrix U defined by 
U = diag{w, it} with u = (n+,n_), we can diagonalize 
Mj_ as 



U^M\ U = 



mo - m 2 ( - E± 
— LAkov 



—\Ako z 
-mo + to 2 ( - 



E\ 



(22) 



Now we can easily show that 

det{Mj_} = [E 2 ± + {Ak) 2 - (m - m 2 C) 2 ] 2 = 0. (23) 

This indicates that there exist doubly degenerate four 
different values of n. Two of which are evanescent solu- 
tions (i.e., k > 0) exponentially damped deep inside the 
sample, while the other two correspond to exponentially 
increasing solutions. By superposing two evanescent so- 
lutions with k = k±, we construct a general solution 
localized near the surface as 



3 K + (X 3 -X sf ) 



U+) 



(24) 



The boundary condition \i/j(x 3 
= \u-) for k + 7^ K-. 



when 



= x s f)) = holds only 
As shown in Appendix 
A, we see that this results in the zero-energy condition 
E± = and 



K±(x 1 , X 2 



A ± y/A 2 + 4TO m 2 
-2m 2 



(25) 
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with A = A + 1TL2A. We also find that two basis eigen- 
states, |+) and |— ), for H± with E± — are given by 



\±)(x\x 2 ,x s ) = p{x 3 ;x\x 2 )\±)){x\x 2 ) 



with 



\±)){x\x 2 ) = -j= 



n+ 
=Fin_ 



and 



p(x 3 ;x\x 2 ) = -±= ( e --(- 3 -^) _ e - + (- 3 -^)) 



(26) 



(27) 



(28) 



where c(x 1 , a; 2 ) is a normalization constant. The x 1 - and 
^-dependences of p arise from c and k±. 

In the third step we derive the effective surface Hamil- 
tonian within the k ■ p-approximation. The following 
derivation is based on the observation that any surface 
state I "J/) can be represented as a linear combination of 
|+) and |— ) with the amplitude respectively specified by 
<5+ and <5_, i.e., 



I*) 



(29) 



Within the k ■ p-approximation, the effective surface 
Hamiltonian H c s for the two-component spinor a — 
*(a+, a_) is given by 



Hefi 



-\H\\\+) (+\H\\\-) 



Here, each matrix element is expressed by 



(a\H {l \a' 



dx 3 VGp((a\H ll \<j'))p, 



(30) 



(31) 



x s f — L c 



where l c is the cutoff length being much longer than the 
penetration depth A. Note that the factor \[Q reflects 
the fact that dV = \ /r Qdx 1 dx 2 dx 3 . Accordingly, we set 
the normalization constant as 

c(x 1 ,x 2 ) = [ ' dx 3 VG (e K - ix3 - Xs!) -e K+{x3 - x " !) ). 

(32) 

Hereafter we use the shorthand notation of the average 
over x 3 defined by 



P Bf da; 3 • • • o 2 

Jx s f — l c >__ 



(33) 



The average in eq. (fH?)) should be identified with this. 

Without loss of generality we hereafter require that n± 
are connected by time-reversal operation as 



(34) 



indicating that if n_ = '(711,-712), then n + = 
(n2,ni). If e 3 is parameterized in terms of spher- 
ical coordinates ^(a; 1 ,^ 2 ) and (^(x 1 ,^ 2 ) as e 3 = 
(sin 9 cos </>, sin 9 sin 0, cos 9) and hence 



cos 9 sin#e 
sin 9e 1( ^ — cos ( 



a natural solution of eq. ([2T]) is given by 
n + (x 1 , x 2 ) 

n_(x 1 ,x 2 ) — 

■ cos 2 



(35) 

(36) 
(37) 



This standard expression of n± obviously satisfies 
eq. (jM)) . Note that in the derivation of the effective sur- 
face Hamiltonian, we require only eq. (I34[) and do not 
explicitly use eqs. (|3"6")l and (|3"T|). 

Let us evaluate the matrix elements given in eq. (|31[), 
We easily find that the diagonal elements vanish, i.e., 
((±|_ff|||±)) = 0. In evaluating the off-diagonal elements, 
we should note that di in H\\ acts on not only |±))(a: 1 , x 2 ) 
but also p(x 3 ; x 1 , a; 2 ). After straightforward calculations 
(see Appendix B), we find 



T-LeS = 



V+ 
V- 



(38) 



with 

2 

^+ = E 

i=l 



(ViA ~ &ma) + - [di (rjiA - &m 2 )] 



(39) 



2 r 



(rjiA - ^m 2 )* di - - [d t faA - &m 2 )*] 



(40) 



where 



dx 3 Vgp 2 n\a l n- 



Xsl—lc 



(VGr 



(VG) 



(41) 



.7 = 1 : •' 



2 ^(,Vgg lJ n\d j n-) 



(42) 



We show in Appendix C that ^ is simplified to 
2 (V^6/n +CT J «-) 



E 



(43) 
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where 60 



(djes) ■ e 1 . The second term in T>± is 



essential in ensuring the hermiticity of T-L c s when the co- 
efficient rjiA — £i?7i2 depends on x' i 17 i 23 i 24 It should be 
noted that the effective velocity in the a^-direction is de- 
termined by rjiA — ^77i 2 , where the second term with ni2 
represents the renormalization due to curvaturei 17 i 20 

In the final step we slightly modify the obtained ef- 
fective Hamiltonian to adapt our framework to the stan- 
dard convention. Since we insert the factor yy in the 
definition of the matrix elements [see eq (|31[) ]. the inte- 
gral measure for the orthonormalization of a. is da; 1 da; 2 . 
However, d(5) = (VG)dx 1 dx 2 is a more natural choice. 
Accordingly we define the new two-component spinor a 
as 

a = , 1 n (44) 



0) 

for which d(S) can be applied. The effective Hamiltonian 
for a is obtained as 



V+ 
£>_ 



(45) 



where 

2 

i=l 
2 



(r H A-Z i m2)(d i + ~[dMVG) 



Ul,A -Lm-2)" [<h~l \(k\n(VQ) 



- [di(rnA-^m 2 )*] 



(46) 



(47) 



The above set of equations, eq. ((35]) with eqs. (|46l) and 
(|4T|) . constitutes the central result of the paper. 

We have thus seen that the obtained effective Hamil- 
tonian takes indeed a generalized Dirac form with linear 
differential operators. Low-energy electrons on the arbi- 
trary curved surface of a topological insulator do obey 
a Dirac equation. It should be also emphasized that 
the obtained linear differential operator involves a cor- 
rection term, (l/2)9Jn(vy), that corresponds precisely 
to what is known as the spin connection in the Dirac 
theory on curved surfaces ^2, ~ 12 i 15 i 25 By noting d(S) = 
(\/G)dx 1 dx 2 , we have identified the precise origin of the 
spin connection as arising from the spatial variations of 
an infinitesimal area element. Last but not the least, we 
have pointed out that the last term in the expression for 
T>± [see eqs. (|46|) and (|47p ] appears for ensuring the her- 
miticity of "Hcff, and is unrelated to the spin connection 
itself. 



IV. ORIGIN OF THE BERRY'S PHASE n 

Let us discuss the origin of the Berry's Phase 7r in 
the context of the boundary condition for a (a; 1 , a; 2 ) = 



*(aq_, a_), and its relation to the so-called spin-to-surface 
locking. Let us consider a situation in which either one or 
both of our curvilinear coordinates x % are cyclic. Angular 
coordinates of a closed (e.g., cylindrical or spherical) ge- 
ometry could be a typical example of such a coordinate. 
In the following we assume that only the coordinate x l 
(i = 1 or 2) is cyclic with a cycle of Li (the other coor- 
dinate does not appear explicitly in the discussion). Our 
starting point is the fact that any wave function 

\*(x i )) = a + (x i )\+)(x i )+a-(x i )\-)(x i ) (48) 

must satisfy the periodic boundary condition, i.e., 

|*(^)) = Wx' + Li)). (49) 

The local spin quantization axis n±(x l ) plays a crucial 
role in our argument. Obviously, if n±(x l ) = n±(x' 1 + Li) 
and hence \±)(x' 1 ) — \±)(x l +Li), the boundary condition 
for a^a; 1 ) must be periodic as 



a(x 



Li). 



(50) 



Note that n± could change its sign as n±(x l ) = 
—n±{x l + Li) after the coordinate x % finishes one com- 
plete cycle of evolution (i.e., x % — > x l + Li). If this is 
the case, the sign of \±}(x l ) is also reversed as \±)(x l ) — 
— \±}(x l + Li). Accordingly, the boundary condition for 
a must be antiperiodic as 



a(x l ) — —a(x l + Li). 



(51) 



This indicates that the boundary condition for a is sim- 
ply determined by whether or not n±(x l ) changes its sign 
when the cyclic coordinate x l is shifted by one complete 
cycle Li (in geometrical terms, this would correspond to 
one complete revolution of the Dirac electron around the 
closed surface). Actually, n±(x l ) changes its sign when it 
rotates by ±2ir around an arbitrary axis in the spin space 
as the Dirac electron revolves once around the closed sur- 
face. 

Let us observe that the antiperiodicity of the bound- 
ary condition discussed above is equivalent to the Berry's 
phase 7r in the Dirac theory on curved surfaces. In the 
latter point of view the antiperiodic boundary condition 
is abandoned (i.e., replaced with the periodic one) at the 
cost of introducing a Berry's phase tt. This can be seen 
as follows: we focus on the case in which n± changes 
its sign as n±(x l ) — —n±(x % + Li). Then, we reformu- 
late this problem by the use of the following single- valued 
basis vectors 



(52) 



h±(xi) = exp ( vk— ) n±(xi), 



which obviously satisfy the periodic boundary condition: 
n±(x l ) = h±(x l + Li). Reflecting the fact that eq. 
dose not hold for n±, the effective Hamiltonian in this 
single- valued basis becomes 



u sv — 

rt cff — 



V% 



(53) 
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where 



n v = E 



(ijiA - &m 2 ) ( di + - diki{VG) 



u 



[di (j]iA - &m 2 )] 



- (rjiA - ^m 2 )* ( di + - di\n(VG) 



- [di (riiA - £im 2 y] 



(54) 



Li 
(55) 



Note that in the above formulas the change of the bound- 
ary condition has been absorbed as a correction to deriva- 
tives, i.e., the term of the form of m / 'Li, which sums up 
to a Berry's phase it. In this sense, the Berry's phase 
7r is a mere rewriting of the antiperiodicity of the basis 
vectors n±, while it is often regarded as an important 
part of the spin connection in literaturesj n i 12 i 16 ' 18 

We have so far argued that the Berry's phase it should 
be attributed to the sign change of the local spin quan- 
tization axis n± caused by a ±2ir rotation in the spin 
space. Previously, the Berry's phase it is interpreted as 
a consequence of the spin-to-surface locking^ -14 ' 16 This 
statement is plausible, but fails to capture the precise ori- 
gin of this term. Indeed, even in the situation where the 
Berry's phase 7r or equivalently the antiperiodic bound- 
ary condition plays a role, the spin-to-surface locking 
does not necessarily or globally occur, as is demonstrated 
in the spherical system of a topological insulator*^ The 
Berry's phase n is encoded in the surface Dirac the- 
ory, specified by an explicit form of the effective surface 
Hamiltonian and by the periodicity of the boundary con- 
dition of n± as well. Contrastingly, the spin-to-surface 
locking is not an ingredient of the surface Dirac theory 
but a characteristic feature of the typical surface wave 
functions, and in this sense not appropriate for being 
given a full status as the origin of the Berry's phase it. 



V. APPLICATION TO REPRESENTATIVE 
CASES 

In this section we apply the general framework estab- 
lished so far to the following two representative cases: 
samples of either a cylindrical or a spherical shape and 
find an explicit formula [corresponding to eqs. (|46[) and 
(|47|) ] of the differential operators T>± that specify the 
explicit form of the Dirac Hamiltonian (|4"5l) . In the anal- 
ysis given below we employ orthogonal curvilinear co- 
ordinates, for which = diag(<?n, (722, 1) and g lJ = 
diag( 5 n , ff 22 ,l) with g 11 = {gu)- 1 and g 22 = (ffa)- 1 . 



A. The cylindrical case 

Let us consider an infinitely long cylindrical topolog- 
ical insulator aligned along the z-axis with radius R. 



We employ the following three coordinates (a; 1 , a; 2 , x 3 ) = 
(4> , z, r), in terms of which the 3D Cartesian coordinates 
are expressed as (x,y,z) — (r cos (j>, r sin^, z). The pa- 
rameter Xgf is simply equal to R. The tangent and nor- 
mal vectors are 



ei = (— rsin</>, r cos0, 0), (56) 
e 2 = (0,0,1), (57) 
e 3 = (cos 0, sin 0,0), (58) 



and 



cos < 



r r 
(0,0,1), 
(cos0, sin0, 0). 



(59) 

e- = (li. I). J i. (60) 

(61) 

The elements of the metric tensors are gn = r 2 and 
g-22 = 1, which results in Q = r 2 , and the coefficients 
of the Weingarten equation are b\ l = — r^ 1 and b 2 — 
&2 1 = b 2 2 = 0. From the expressions of e l we obtain the 
spin matrices as 





1 
-1 

e~ u 




(62) 
(63) 
(64) 



As the unit vectors satisfying a 3 n± = ±n±, it is conve- 
nient to use those given in eqs. (|36|) and (|37|) at 9 = ir/2, 



n± 



1 

71 



e 2 

■ 4> 

ie 1 * 



(65) 



Then we immediately find that n+^ri- = ir 1 and 
nt_cr 2 n_ = 1. Substitution of these results with \[Q = r 



and 61 1 = —r 1 and b\ 2 
and 103]) yields 



b, 1 
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(VG) 



= 1, 



-(- 

(r) \ r 



= 0. 



into eqs. ((41 



(66) 
(67) 
(68) 
(69) 



Noting that e^ln(vy) = d z ln(^/Q) = we finally obtain 
the differential operator T>± as 



P±=iU 



m 2 



(r) 



±Ad z 



(70) 
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If the penetration depth A for surface states is much 
shorter than R, we can approximate as (r) = R and 
(r _1 ) = Rr 1 , and then T>± is simplified to 



A 



(71) 



The effective Hamiltonian is given by eq. (j4"5"j) with I?± 
obtained above. The result similar to this has been re- 
ported in refji£, where the renormalization correction 
m 2 jR to the effective velocity is ignored. Let us con- 
sider the boundary condition for a spinor wave function 
a(<j),z). Since n±(<p) in eq. (|65|) changes its sign when 



cj) — > <p + 27r, the boundary condition for the variable 
must be antipcriodic, i.e., 



< + 2w,z). 



(72) 



B. The spherical case 

We turn to the second case of a spherical topologi- 
cal insulator with radius R. We employ the standard 
spherical coordinates (a; 1 ,^ 2 ,^ 3 ) = (9,(f),r), in terms 
of which the 3D Cartesian coordinates are expressed as 
(x, y, z) = (r sin # cos rsin#sin0, rcos9). The param- 
eter x a { is again equal to R. The tangent and normal 
vectors are 



ei = (r cos 9 cos 0, r cos 6 sin <fi, —r sin ( 
e 2 = (— r sin # sin 0, r sin cos 0, 0), 
e3 = (sin 9 cos (j>, sin sin 0, cos 9) , 



and 



e = 



cos 9 cos ^ cos 9 sin < 



sin6> 



e = - 



r 

sin c 



e = 



r sin t/ r sin f 
(sin 9 cos (/), sin 9 sin 0, cos ( 



(73) 
(74) 
(75) 



(76) 

(77) 
(78) 



that nta^re 



_1 and ret er 2 n_ 



Sub- 



stitution of these results with = r 2 sin#, &1 1 = b 2 2 = 



and6i 2 



6 



= &2 1 = into cqs. (gT|) and (g3j yields 

(r- 1 ^) (r 



i((rsin0)- 1 \/g) 



r- 2 VG) _ 

(VQ) 

-i((r 2 sin# 



(r 2 ) sin# 



1 



(VG) 



(82) 
(83) 
(84) 
(85) 



Noting that <%ln(V£) = cot 6 and <9 ln(\/Cj) = we fi- 
nally obtain the differential operator T>± as 



d,<, 



l 

=F - cot I 
2 



(86) 



It should be emphasized that (1/2) cot# is identified with 
the spin connection in the curved Dirac theory*^* 2 - 5 If the 
penetration depth A for surface states is much shorter 
than R, we can approximate as (r) = R and (r 2 ) = R 2 , 
and T>± is simplified to 



A 



m 2 \ 1 
~RV R 



0,i 



T \ cot 9 



(87) 



The effective Hamiltonian is given by eq. (|45|) with V± 
obtained above. The result similar to this has been re- 
ported in reft^, where the renormalization correction 
m 2 /R to the effective velocity is ignored. Let us con- 
sider the boundary condition for a spinor wave function 
a(9,(j)). The system is periodic with respect to </>, so we 
consider the boundary condition when <p is increased by 
2ir. Since n±(cfi, 9) changes its sign when </> — > <P + 2tt, the 
boundary condition for the variable cj) must be antiperi- 
odic, i.e., 



-a(9, cj) + 2?r) 



The elements of the metric tensors are gn — r 2 and 
g 2 2 = t 2 sin 2 9, which results in Q = r 4 sin 2 9, and the 
coefficients of the Weingarten equation are foi 1 = b 2 2 = 
— r _1 , and bi 2 = b 2 x = 0. From the expressions of e 1 we 
obtain the spin matrices as 



r 



— sin 9 cos 9e 

cos 6>e'^ sin 9 

-e-^ 

e 1 * 



rsin# 

cos 9 sin 9c" 
sin 9e lc ^ — cos ( 



(79) 
(80) 
(81) 



As the unit vectors satisfying <r 3 n± = ±n±, it is conve- 
nient to use those given in eqs. (|3l)]) and (|3"7|) . We find 



VI. SUMMARY AND DISCUSSION 

The behavior of low-energy electrons on an arbitrary 
curved surface of 3D (strong) topological insulators has 
been considered on general grounds. In contrast to the 
specific cases studied earlier, we have reached a unified 
description of such low-energy electrons by giving the 
most general form of the surface Dirac Hamiltonian [eq. 
(|4"5j) with eqs. fH]) and (|37]1] that has been explicitly 
derived from the bulk effective theory in the continuum 
limit. It was shown that the low-energy surface electrons 
do obey the Dirac equation in this generalized form with 
the effective velocity renormalized by the curved nature 
of the surface. A special attention has been paid to the 
boundary condition for a spinor wave function a, which 
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becomes relevant on a closed surface described by at least 
one cyclic coordinate x 1 . Whether the boundary condi- 
tion is periodic or antiperiodic depends on the behavior 
of the local spin quantization axis n±. If the sign of n± 
is unchanged after one complete cyclic evolution of the 
coordinate x l , the boundary condition for a is periodic. 
On contrary, if n± changes its sign due to a ±2n rota- 
tion in the spin space, the boundary condition becomes 
antiperiodic. It is argued that the antiperiodicity of the 
boundary condition is equivalent to the Berry's phase tt. 

Previously, the effective Hamiltonian for Dirac elec- 
trons on a curved surface has been considered in a 
framework different from the one presented in this pa- 
per ^ r 12 ' 15 In this alternative viewpoint, one starts from 
a two-dimensional Dirac equation for a fiat surface, and 
takes account of the curved nature of a surface by a co- 
ordinate transformation, resulting in the curved surface 
Dirac theory. The effective Hamiltonian thus obtained 
contains a fictitious vector potential called the spin con- 
nection, which corresponds to the term (l/2)djln(\/5) in 
our framework. We have shown that the spin connec- 
tion represents corrections arising from the spatial varia- 
tion of an infinitesimal area element. We have also seen 
that such an ad hoc curved Dirac theory overlooks the 
renormalization of the effective velocity arising from the 
quadratic mass term 777,2 • This is because the theory ig- 
nores from the outset the three dimensional nature of the 
problem. It should be noted that the above renormaliza- 
tion arises even in the limit in which the penetration 
depth A of surface states is vanishingly short, as is seen 
in 85. 
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Appendix A: Proof of E± 







From the above equation we easily observe that if 
the following two equations 



mo - in 2 (+ - Ej_ too - m 2 (- - E± 



Ak+ Ak- 
m-o - m 2 C+ + Ex too - TO 2 C- + E± 



An, 



Ak_ 



(A3) 
(A4) 



must hold simultaneously. This directly results in E± = 
0. Under this zero energy condition, det{M^} = yields 
to — to 2 C± = Ak±, where to > and 7712 < are 
assumed. Solving mo — TO2C± = Ak± with respect to k±, 
we obtain eq. ([25]) . 



Appendix B: Derivation of the Off-Diagonal Matrix 
Elements 



From cqs. fTTjl. @ and (J27J) it is easy to show that 



\H n 



X„f — lc 



dx 3 Vgp[Aj2 

n' + a' l d i 'n_ 



— TO 2 «.^A||n_ J p. (Bl) 

Let us denote the first and second terms in the right- 
hand side of eq. (|B1I) as )i and (+|ifii|— ) 2 , re- 
spectively. The first term is rewritten as 



(+\H\\ 



E 



i=l 

where rji is defined in eq. (|41l) and 

1 r x »t 



mAd, + - {d lVl A) + APi 



(B2) 



Pi = -: 



dx 3 p 2 



X B S — lc 



(B3) 



We can show (din} + )a l "n- — n\_a l {dinJ) — by using 
eq. (|3~4"|) . and Pi is reduced to 



As noted in the text, the boundary condition \ip(x 3 = 
x s {)) = holds only when \u + ) = |7t_) for n + 7^ Here 
\u±) are the eigenvectors satisfying M±(k±)\u±) = 0, 
where 



Mj_(k±) 



to 



to 2 C± - E± 
-L4k±ct 3 



-m 



-iAn±a 3 
f to 2 C± - 



(Al) 



with C± — K ± + A«±. It is instructive to rewrite the 
eigenvalue equation as 



nip — m 2 Q±-E± 
Ak± 

-i<7 3 



— 1(7 

mo-m 2 (±+E ± 
Ak± 



\u±) = 0. 



(A2) 



Pi = - : . 



X s f—l c 



dx 3 p 2 n t + (d i VG(T i )n. 



(B4) 



It is easy to show with e 1 = e 2 X e^/^/Q and e 2 = e% x 
ei/VG that 

dx^/Ga 1 = [(ftea) x e 3 + e 2 x (ftea)] • <x, (B5) 
d 2 V<3o- 2 = [(d 2 e 3 ) x ei + e 3 x (ftei)] • a. (B6) 

Applying the Gauss equation (TTTT) and the Weingarten 
equation (fT4|) . we find that 



dxVGcr 1 =Vg(T% 1 a 1 -Tl 1 a 2 + b l 1 o- 3 ), (B7) 
d 2 VQ(J 2 = VG {-T^a 1 + T\ 2 o 2 + 6 2 V) . (B8) 
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Substituting these into eq. (|B4[) and using = and 

n+cr 3 n_ = 0, we find that Y^h=i Pi — 0- Now we turn to 
the second term which is given by 

2 



|ff|||->2 = -m 2 J2 I ' d *W 



d l (vW^) n^p. (B9) 



Since n+n_ = 0, only the terms with n^_<9jn_ or 
nldidjU- do not vanish. It is then reduced to 



\Hn\-h = Yl 



-S,im 2 di - - (di^im 2 ) + m 2 Qi 



(BIO) 



where & is defined in eq. (|42l) and 



Combining the resulting (+|JJy|— )i and ) 2 we 

finally arrive at eq. (f59"|). The expression of (— |ifii |+) can 
also be obtained by repeating the procedure described 
above. 



Appendix C: Simplification of fi 

In this short Appendix we simplify the expression of 
£i defined in eq. (|42|) . The starting point is the following 
eigenvalue equation: er 3 n_ = — n_. Differentiating this 
by x 3 and then constructing the inner product between 
the resulting expression and n + , we obtain 



nl(djnJ) = -~n\. (<V 3 )n_ 



(CI) 



The Weingarten equation ([M)) enables us to replace 



dj<j 3 = (dje 3 ) ■ cr with — J2k=i bjk& ■ This results 



.7 = 1 



d^V^pV^^riD^n-). (Bll) 



6=5: 



da; 3 VSp 2 g ij b 3 ; fc n\ a k n 



(C2) 



We can show with eq. (IMf that 

E?.,-=i^(* n +)( fi, i n -) = °- Hence TLiQi = °- 



Noting that Ylj=i 9 io bjk — bk , the expression (|C2[) for 
is reduced to eq. p3)) . 



1 L. Fu, C. L. Kane, and E. J. Mele: Phys. Rev. Lett. 98 
(2007) 106803. 

2 J. E. Moore and L. Balents: Phys. Rev. B 75 (2007) 
121306. 

3 R. Roy: Phys. Rev. B 79 (2009) 195322. 

4 D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, 
and M. Z. Hasan: Nature 452 (2008) 970. 

5 Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. 
Bansil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. Hasan: 
Nat. Phys. 5 (2009) 398. 

6 Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. Mo, 
X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. C. 
Zhang, I. R. Fisher, Z. Hussain, and Z.-X. Shen: Science 
325 (2009) 178. 

7 A. Nishide, A. A. Taskin, Y. Takeichi, T. Okuda, A. Kak- 
izaki, T. Hirahara, K. Nakatsuji, F. Komori, Y. Ando, and 
I. Matsuda: Phys. Rev. B 81, 041309 (2010). 

8 K. Kuroda, M. Ye, A. Kimura, S. V. Eremeev, E. E. 
Krasovskii, E. V. Chulkov, Y. Ueda, K. Miyamoto, T. 
Okuda, K. Shimada, H. Namatame, and M. Taniguchi: 
Phys. Rev. Lett. 105 (2010) 146801. 

9 M. Z. Hasan and C. L. Kane: Rev. Mod. Phys. 82 (2010) 
3045. 

10 D.-H. Lee: Phys. Rev. Lett. 103 (2009) 196804. 

11 Y. Zhang, Y. Ran, and A.Vishwanath: Phys. Rev. 79 
(2009) 245331. 

12 Y. Zhang and A.Vishwanath: Phys. Rev. Lett. 105 (2010) 
206601. 



1 R. Egger, A. Zazunov, and A. Levy Yeyati: Phys. Rev. 
Lett. 105 (2010) 136403. 

14 J. H. Bardarson, P. W. Brouwer, and J. E. Moore: Phys. 
Rev. Lett. 105 (2010) 156803. 

15 V. Parente, P. Lucignano, P. Vitale, A. Tagliacozzo, and 
F. Guinea: Phys. Rev. B 83 (2011) 075424. 

16 K.-I. Imura, Y. Takane, and A. Tanaka: Phys. Rev. B 84 
(2011) 195406. 

17 Y. Takane and K.-I. Imura: J. Phys. Soc. Jpn. 81 (2012) 
093705. 

18 K.-I. Imura, Y. Yoshimura, Y. Takane, and T. Fukui: Phys. 
Rev. B 86 (2012) 235119. 

19 K.-I. Imura, M. Okamoto, Y. Yoshimura, Y. Takane, and 
T. Ohtsuki: Phys. Rev. B 86 (2012) 245436 

20 K.-I. Imura and Y. Takane: larXiv:1211.2088l 

21 C.-X. Liu, X.-L. Qi, H. Zhang, X. Dai, Z. Fang, and S.-C. 
Zhang: Phys. Rev. B 82 (2010) 045122. 

22 W.-Y. Shan, H.-Z. Lu, and S.-Q. Shen: New J. Phys. 12 

(2010) 043048. 

23 A. Raoux, M. Polini, R. Asgari, A. R. Hamilton, R. Fazio, 
and A. H. MacDonald: Phys. Rev. B 81 (2010) 073407. 

24 R. Takahashi and S. Murakami: Phys. Rev. Lett. 107 

(2011) 166805. 

25 A. A. Abrikosov: Int. J. Mod. Phys. A 17 (2002) 885. 



